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1. Introduction 



In the traditional approach to quantum field theory, one starts from classical fields and 



a Lagrangean which describes the interaction. These objects get quantized and S'-matrix 
elements or Greens functions are constructed with the help of the Feynman rules. A typical 
example which is used in this paper as a model theory is the massless <3? 3 -theory, where 
the interaction Hamiltonian density is given by the normal ordered third order monomial 
• of a free uncharged massless scalar field and a coupling constant A 

§ n mt (x) = A : $ (x )3 (1) 



and x is an element of Minkowski space, denoted by K, 4 in the following. The perturbative 
(]} \ S-Matrix is then constructed according to the expansion 

*^ - ( - / )" /" 

S = l + Y, / dx\...dxiT{H in t(x x )H in t(x 2 ) • ... • n int (x n )}, (2) 

where T is the time-ordering operator. It must be pointed out that the perturbation 
series eq. (J2J) is formal and it is difficult to make any statement about the convergence 
of this series. Furthermore, two problems arise in the expansion given above. First, the 
time-ordered products 

T n (xi,X 2 , ...,X n ) = {-l) n T{Hint{xi)Hint{x2) ■ ■■■ ■ Hint{x n )} (3) 

are usually plagued by ultraviolet divergences. However, these divergences can be removed 
by regularization, such that the operator-valued distributions T n can be viewed as well- 
defined, already regularized expressions. Second, infrared divergences are also present in eq. 
(J2J- This is not astonishing, since the T n 's are operator- valued distributions, and therefore 
must be smeared out by test functions in 5(IR, 4n ), the Schwartz space of functions of rapid 
decrease. One may therefore introduce a test function g(x) € 5(]R, 4 ) which plays the role 
of an 'adiabatic switching' and provides a cutoff in the long-range part of the interaction, 
which can be considered as a natural infrared regulator Then the infrared regularized 

S'-matrix is given by 

00 1 f 

S(g) = 1 + ~ / dxf...dx n T n (xi, x n )g(x l ) • ... • g(x n ), (4) 

n=l n ' 
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and an appropriate adiabatic limit g — > 1 must be performed at the end of actual calcu- 
lations in the right quantities (like cross sections) where this limit exists. This is not one 
of the standard strategies usually found in the literature, however, it is the most natural 
one in view of the mathematical framework used in perturbative quantum field theory. 
The vertex diagram appears in $ 3 -theory at third order of perturbation theory via 

T 3 ( Xl ,x 2 ,x 3 ) = (^pj T{: ^( Xl ) :: <S> 3 (x 2 ) :: <D 3 (x 3 ) :} 

= A 3 V(3;i,X2,X3) : ^(xi)^(x 2 )^(x 3 ) : +other graphs, (5) 

where V{x\, X2, x 3 ) = Ap(x± — x 3 )Ap(x 3 — X2)Ap(x2 — x\) is given by the scalar Feynman 
propagator discussed below in further detail. Of course, the massless vertex functions plays 
also a prominent role in applied quantum field theories like QCD, where it emerges in the 
three-gluon vertex or in perturbative models of quantum gravity. 

In this paper, the Fourier transform of the massless vertex distribution (also called 
vertex function in the sequel) 

V(p,q)= [ d 4 x 1 d 4 X2V(x 1 ,X2,x 3 )e ip( ^ Xl - x ^ +i ^ x;i - x ^ 



4 I [(p~k) 2 + iO] [(q - k) 2 + iO] [k 2 + iO] (6) 

will be calculated in a non-standard manner by making explicit use of causality as a 
fundamental property of local quantum field theories 3|. The presented dispersive method 
is based on considerations in configuration space and integrals in momentum space are 
not evaluated by the help of Wick rotations. Furthermore, no mass terms are needed as 
regulators of infrared divergences, the vertex function is constructed directly without any 
limiting procedure. 

It should be noted that the vertex function is translation invariant in the sense that 
V(xi,X2,x 3 ) = V(x\ — a, X2 — a, x 3 — a), where a is an arbitrary four- vector. This symmetry 
property can be inferred directly from the decomposition of the vertex function into three 
Feynman propagators given above. Introducing new integration variables z\ = x\ — x 3 and 
Z2 = x 3 — X2, the Fourier transform in eq. (jHJ) can also be written as (Ap(x) = Ap(—x)) 

V(p, q)=J d i z 1 d i z 2 V{z x + S3, ~Z2 + x 3 , x 3 )e^ +l ^ 

= J d 4 Zl d 4 Z2 A F (z 1 )A F (z 2 )A F (z 1 + z 2 y vz ^ iqz \ (7) 

which shows that V(p, q) does indeed not depend on x 3 . 

The focus of this work is on the causal and analytic properties of the vertex function and 
its infrared-finite dispersive calculation in the massless case. The more involved general- 
mass case of the vertex function has been treated by 't Hooft and Veltman jS], such that 
results for the massless case could be obtained by considering the appropriate limit. An 
explicit result for the massless case has been given in and a discussion of the vertex 
function within a special framework based on a geometrical interpretation of kinematic 
invariants and other quantities of one-loop diagrams is contained in [7j, where a short 
overview over different calculational methods for one-loop N-point functions can also be 
found. 
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Before investigating the causal structure of the scalar vertex function, some important 
analytic properties of the Feynman propagator and other related tempered distributions 
will be recalled in the following section. 



2. Properties of tempered distributions 

The Feynman propagator of a free scalar and charge neutral quantum field $(x) fulfilling 
the wave equation 

□ $( x ) = d^^ix) = (8) 
is given in configuration space by the vacuum expectation value 

A F (x) = A F ( Xl - x 2 ) = -»<0|T(#(a?)#(0))|0) = -s(0|T($(xi)$(x 2 ))|0) 

d ' k e ~ tkx = A L_ = ^pl _ ^ 5(x2) (9) 

(2tt) 4 fc 2 +i0 4vr 2 x 2 -i0 4vr 2 x 2 4ir K h K > 

where P denotes principal value regularization and 5 is the one-dimensional Dirac distri- 
bution depending on x 2 = x^x^ = (x ) 2 — (x 1 ) 2 — (x 2 ) 2 — (x 3 ) 2 = x\ — x 2 (Hj. 

Knowledge of the basic properties of the Feynman propagator is sufficient for many 
applications in perturbative quantum field theory. For the forthcoming discussion, it is in- 
structive to review some of the useful properties of basic tempered distributions appearing 
in quantum field theory. A scalar neutral field can be decomposed into a negative and 
positive frequency part according to 

$(x) = $-( s ) + = / __^ = a(k)e~ ikx + a\k)e ikx (10) 

J V(2vr)3 2 |fc| 

with 

[a(k),a\k')} =5^{k-k'), [a(k),a(k')] = [a\k),a\k')]=0. (11) 

The commutation relations for such a field are given by the positive and negative frequency 
Jordan-Pauli distributions 

A ± {x) = -*[$ T (x), $±(0)] = -i(0p T (x), ^(0)1)0) , (12) 

which have the Fourier transforms 

A ± {k) = J d 4 xA ± (x)e ikx = T(2vri) @(±k°)5(k 2 ). (13) 

The fact that the commutator 

[$(x),$(o)] = iA + (x) +iA~(x) =: iA(x) (14) 
with the Fourier transform 

A(Jfe) = -(2iri)sgn{k°)5{k 2 ) (15) 

vanishes for spacelike arguments x 2 < due to the requirement of microcausality, leads 
to the important property that the Jordan-Pauli distribution A has causal support, i.e. it 
vanishes outside the closed forward and backward light cone such that 

suppA(x) C F~ UV + , V ± = {x\x 2 > 0, ±x° > 0} (16) 
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in the sense of distributions. A further crucial observation is the fact that one can introduce 
the retarded propagator A ret (x) which coincides with A(x) on V — {0}, i.e. (A ret ,(p) = 
(A, ip) holds for all test functions in the Schwartz space (/3G5(]R, 4 ) with support suppt/? C 
K, 4 — V . One may write in configuration space 

A ret {x) = Q{x°)A(x), (17) 

and transform this 'splitting' formula into a convolution in momentum space 1 

^ ret (k)= [ J^P A(p)6(fc-p). (18) 



, (2^) 4 

The Fourier transform of the Heaviside distribution Q(x°) can be calculated easily 

e(ifc) = lim [ d A x@(x )e- tx0 e ikox °- ilis = §^ 5^ Ik). (19) 

e^OJ fe u + i0 

For the special case where k is in the forward light cone V + , one can go to a Lorentz frame 
where k = (k°, 0) such that eq. (fl8|) becomes 

Hence, for arbitrary k £ V + , A ret would be given by the dispersion relation 



A ret tk) = — dt . (21) 

K ' 2ttJ 1-t + iO K J 

However, the integral in eq. (|21|) is undefined in the massless case considered here. This 
illustrates the fact that it is not possible in general to convert the product of two distri- 
butions into a convolution via a Fourier transform, even though expression (|17jl describes 
a well-defined distribution. Here, one can circumvent this problem, e.g., by introducing a 
mass for the field $ as a regulator for the corresponding massive Jordan-Pauli distribution 
which is also defined as a field commutator (for a different method, see, e.g., For the 
massive Jordan-Pauli distribution A m (x) one gets in momentum space 

A m (k) = -(2TTi)sgn(k°)5(k 2 - m 2 ) (22) 

and the following expression for the massive retarded distribution Aj^*(/c) (k£ V + ) 

sgn(tk°)8(t 2 k 2 - m 2 ) 



t + iQ 



2VFm(l-t + iO) (23) 

The ^-distribution in the Fourier transformed Jordan-Pauli distribution eq. (|22|) contains 
a mass term which is not present in eq. (|15jh such that A m (k) does not exhibit a singular 
behavior on the light cone in momentum space. 

1 The Heaviside distribution O(x ) could be replaced by Q(vx) with an arbitrary vector v G V + = 
V + — dV + inside the open forward light cone. 
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As a special case of the edge of the wedge theorem jl] it is known that the Fourier 
transform of the retarded distribution A ret (x) is the boundary value of an analytic function 
r(z), regular in T + := R 4 + iV + ^E3- This way one obtains from r(z) = l/(z 2 — m 2 ), 
z 6 T + , and m —* 

The analytic expression for the Feynman propagator is then recovered, and one obtains 

A re \k)=A F (k) + A-(k). (25) 

It is helpful for the forthcoming to understand eq. I|25|). This is best done in configura- 
tion space. One obviously has 

iA F (x x -x 2 ) = <0|r($(xi)$(s 2 ))|0) 

= (0|G(x? - ^)$(si)$(x 2 )|0) + (0\@(x° 2 - x?)$(x 2 )$(xi)|0) (26) 

and 

iA~(xi - x 2 ) = -(0|$(x 2 )$(xi)|0> 
= -<0|Q(x° - x$)*(x 2 )*(x 1 )|0) - (0|Q(x? - x°)*(x 2 )*(x!)|0), (27) 
and therefore eq. (|24() is fulfilled 

iA re *(xi - x 2 ) = iG(x? - x° 2 )A( Xl - x 2 ) 

= (0|Q(x? - x^)$(xi)$(x 2 )|0) - (0|e(x? - x^)$(x 2 )$(xi)|0). (28) 

The discussion presented here can be applied the same way to second order loop diagrams 
^2 El EH- The problem of ultraviolet divergences can be solved in such cases by using 
subtracted dispersion relations pQ. 



3. Causality 

The causal properties of the Jordan-Pauli distribution and the retarded propagator can be 
generalized in an analogous manner to the case of the vertex function. For this purpose, one 
defines the causal operator-valued three-point functions Af op , A^ p and related C-number 
valued contractions A^ and A3, which are related to the vertex diagram via 

A3 }0 p(xi, X 2,x 3 ) = -T 2 (x 2 ,x 3 )Ti(xi) - Ti(x 3 )T 2 (xi,x 2 ) - T 2 (xi, x 3 )Ti(x 2 ) 

+r 1 (x 3 )r 1 (x 1 )T 1 (x 2 )+T 1 (x3)r 1 (x 2 )r 1 (xi) 

(xi,x 2 ,x 3 ) : <I>(xi)<I>(x 2 ) < I ) (x3) : +other graphs, (29) 
A^ op (xi,x 2 ,x 3 ) = +Ti(xi)T 2 (x 2 ,x 3 ) + T 2 (xi,x 2 )Ti(x 3 ) + T 1 (x 2 )T 2 (x 1 , x 3 ) 
-Ti(xi)Ti(x 2 )Ti(x 3 ) - Ti(x 2 )Ti(xi)Ti(x 3 ) 
= A 3 A3 (xi, x 2 , x 3 ) : $(xi)$(x 2 )$(x 3 ) : +other graphs, (30) 



and 



Ag p (xi,x 2 ,x 3 ) = AT (xi,x 2 ,x 3 ) + A 3 (xi,x 2 ,x 3 



6 



ANDREAS ASTE 



= [T 1 (x 1 ),T 2 (x 2 ,x 3 )} + [T 1 (x 2 ),T 2 {x 1 ,x 3 )} - [Ti(x 3 ), T 2 (x x , x 2 )] 
+ [T 1 (x 3 ),T 1 (x 1 )T 1 (x 2 )] + [T 1 (x 3 ),T 1 (x 2 )r 1 (xi)] 
= A 3 A 3 (xi,x 2 ,x 3 ) : $(xi)$(x 2 )$(x 3 ) : +other graphs, (31) 
where Ag P and A 3 have causal support as distributions in the sense that 

suppA(° p) (x 1 ,x 2 ,x 3 ) CT+UTg, (32) 

where 

T± = {(xi,...,x n )\x n GR^, (xj - x n ) £V Vj = 1, ...,n- 1}. (33) 

Therefore, only when both arguments x\ and x 2 are in the forward light cone or both 
arguments x\ and x 2 are in the backward light cone with respect to x 3 , (xi, x 2 , x 3 ) belongs 

to the support of Ag° p \ For the sake of brevity, only a short outline how one can prove 
the special support property of Ag P is given here by considering the case where, e.g., 
x\ £ V + (x 3 ) and x 2 V + (x 3 ). In this case, one can assume that one is in a Lorentz 
system where x^ > x|] > x®, since the transformation properties of all T n 's are fixed by the 
representation of the Lorentz group on the Fock space of massless scalar particles. Then 
one has, e.g., T 2 (x 2 ,x 3 ) = T(Ti(x 2 )Ti(x 3 )) = Ti(x 3 )Ti(x 2 ), and accordingly for the full 
expression eq. (|3*Tj) 

A 3 p (xi,x 2 ,x 3 ) = 

+Ti(xi)Ti(x 3 )Ti(x 2 ) - Ti(x 3 )Ti(x 2 )Ti(xi) 

+Ti(x 2 )Ti(xi)Ti(x 3 ) - ri(xi)Ti(x 3 )Ti(x 2 ) 

-r 1 (x 3 )T 1 (x 1 )T 1 (x 2 )+T 1 (x 1 )T 1 (x 2 )r 1 (x 3 ) 

+T 1 (x 3 )T 1 (x 1 )T 1 (x 2 ) - T 1 (x 1 )T 1 (x 2 )T 1 (x 3 ) 

+ r 1 (x 3 )T 1 (x 2 )T 1 (x 1 ) - r 1 (x 2 )T 1 (xi)T 1 (x 3 ) = 0. (34) 

The same way one may check for all other cases that Ag P vanishes whenever the causal 
condition (|32j) is not fulfilled. The general definition of causal distributions A° p (xi, x n ) 
at higher orders can be found in . 

One should point out that x 3 plays a distinguished role in the discussion above. Whereas 
the time-ordered distribution T 3 (xi,x 2 ,x 3 ) = T 3 (xi 1 ,Xi 2 ,Xi 3 ) is invariant under a permu- 
tation (1,2,3) — > (ii,i 2 ,i 3 ) of the arguments according to its definition in eqns. 1I2I3I4I) . 
Ag P (xi, x 2 , x 3 ) is not. However, the construction presented above would go through with 
xi or x 2 as distinguished arguments. The crucial point is that the causal support property 
of Ag P and its C-number part A 3 will enable the calculation of a distribution which is 
retarded with respect to x 3 via a double dispersion relation in momentum space, as will 
be shown in the following. This retarded distribution Ag 3 * can be related directly to 
the vertex function in analogy to eq. (|25|) . The vertex function again possesses the full 
permutation symmetry in configuration space by construction. 

The vertex function will be constructed below by performing the following calculational 
steps. The Fourier transform A 3 (p, q) = J d 4 xid 4 x 2 A 3 (xi, x 2 , x 3 )e ip<yXl ~ x ' i ' >+ ' iq ^ Xi ~ X2 \ which 
is independent of x 3 due to translation invariance (see also chapter 3.1 in pP), is calculated 
from the expression eq. ()31|) in configuration space. Then the corresponding retarded dis- 
tribution Ag^xi, x 2 , x 3 ) = 0(xi — x 3 )G(x 2 — x 3 )A 3 (xi, x 2 , x 3 ) is evaluated in momentum 
space by a dispersion relation analogously to eq. (|21|h and finally the vertex function is 
obtained directly from analytic considerations. It is also straightforward to check that 
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on the operator level, the time-ordered product T 3 can be obtained from the retarded 
(appropriately regularized) distribution A3 P ' re *(xi, x 2 , x 3 ) 

T 3 (x 1 ,x 2 ,x 3 ) = A° 3 p ' ret (x 1 ,x 2 ,x 3 ) - A 3 op (x 1 ,x 2 ,x 3 ), (35) 

i.e. for the vertex distribution holds 

V(xi,x 2 ,x 3 ) = A r 3 et (xi,x 2 ,x 3 ) - A3 (x 1 ,x 2 ,x 3 ) (36) 

in close analogy to eq. (|25j). 

4. Calculation of A3 

Starting from the first and second order tree diagram given by 

T 2 ( Xl ,x 2 ) = (-if^T : <& 3 (si) :: $ 3 (x 2 ) := 
A 2 

— i— : & 2 (xi)& 2 (x 2 ) : A F (x\ — x 2 ) + other graphs, (37) 



one readily obtains for A3 in configuration space by exploiting the relevant Wick contrac- 
tions in eq. (|3*TT) 





A 3 (xi 


,x 2 ,x 3 ) = Aj(x x , 


x 2 ,x 3 ) + A 3 


(xi,x 2 ,x 3 ) -- 






+ A + (X! - 


-x 2 )A + ( Xl 


- x 3 )A F {x 2 - x 3 ) 


-A+(x 2 - 


xi)A+(x 3 - 


xi)A jF (x 2 


- x 3 ) 


+A+(x 2 


- xi)A + (x 2 


- x 3 )A F (x 1 - x 3 


) -A+(xi- 


x 2 )A+(x 3 - 


x 2 )A F (xi 


- x 3 ) 


+A F (xi 


- x 2 )A + (xi 


- x 3 )A + (x 2 - x 3 


) -A F ( Xl - 


x 2 )A+(x 3 - 


xi)A+(x 3 


- x 2 ) 


+A+(x 3 


- xi)A + (x 3 


- x 2 )A + {xi - x 2 


) -A+(a?i- 


x 3 )A+(x 2 - 


x 3 )A+(xi 


- x 2 ) 


+A+(x 3 


- x 1 )A + (x 3 


- x 2 )A + (x 2 - x\ 


) -A+(xi- 


x 3 )A+(x 2 - 


x 3 )A+(x 2 


-xi).(38) 



Note that a Wick contraction without time-ordering leads to a A + -distribution instead of 
a Feynman propagator. Making use of A + (— x) = — A~(x) and introducing the advanced 
distribution A av such that A(x) = A ret (x) - A av (x) and A ret (-x) = -A av (x), eq. ((3SJ) 
can be simplified to 

A 3 (xi,x 2 ,x 3 ) = 

+A~(xi - x 3 )A+(x 3 - x 2 )A F (xi - x 2 ) + A + (xi - x 3 )A r ' e *(x 3 - x 2 )A+(xi - x 2 ) 
+ A - (xi - x 3 )A av (x 3 - x 2 )A + (x 2 - xi) - {(xi,x 2 ,x 3 ) <-» (-x 2 , -xi, -x 3 )}. (39) 
Performing the Fourier transform, one obtains 

A 3 (p,g) = / d 4 x 1 d 4 x 2 A(x 1 ,x 2 ,x 3 )e ip(xi - X3)+iq{x3 - X2) = (2vr)~ 4 / d*kx 

[A-(p-k)A + (q-k)A F (k)+ 
A+(p - k)A ret {q - k)A + {k)+ 
A~(p - k)A av {q - k)A + {-k)} -{p^ q}. (40) 



8 



ANDREAS ASTE 



From the integral eq. (|4U|) one immediately derives the following symmetry properties of 

A 3 (p, q) = -A 3 (q,p) = A 3 (-p, -q), 

A 3 (-p,q)=A 3 (-q,p). (41) 
As an illustration, the full calculation of the first Fourier integral appearing in eq. (|4Uj) 

(2ir)- 2 h(p,q) = (2^)- 4 | d 4 k[A~ (p - k)A + (q - k)A F (k)} 

= (2n)- 2 J d*kQ(k -p°)5((p - k) 2 )e(q° - k°)5((q - k) 2 )^-^ (42) 

is presented here. Introducing P = p — q and using a new integration variable k' = k — p, 
one obtains 

h(p,q) = J {k , J* - iQ Q(k' )5(k' 2 )e(-P° - k'°)S((P + k'f). (43) 

The first G and 5 distributions can be written as S(k° — \k\)/2\k\ (k' is replaced by k in the 
following). Taking into account the second G distribution one easily derives P° < — k < 0, 
otherwise integral above would vanish. For spacelike P, P 2 < 0, there exists a Lorentz 
frame with P° > such that I\ also vanishes. For timelike P, there exists a Lorentz frame 
where P = = p— q. The last 6 distribution in eq. (|43[) then implies (P°+k°) 2 — k 2 = 0, and 
combined with k° = \k\ from the S distribution 5(k° — \k\) one can replace k° by -P°/2 
in the Feynman propagator appearing in eq. (|43l) after performing the trivial integral 
over kP. For timelike P, there remains to calculate with <5((P° + k ) 2 — k 2 ) replaced by 
5(P 2 + 2P \k\) = S(\k\ -9)/2|P°| 

Hp q) = e(-P°) / ** m - lJ?) (44) 

with k in the denominator above given by k = (-P°/2, k), such that one can write 

(k + p) 2 = (k° + p ) 2 - (k + p) 2 = (-ip° + p ) 2 - P - p 2 - 2|jfe| |p | cos i?, (45) 



and one has |£| = |P°|/2. Hence (P° = P 

+: 

h( P ,q) = Y e(-p°)e(p 2 ) / 



+i 

27r „ , /" dcost? 



J (-P°/2 + p ) 2 -k 2 -p 2 - 2\k\ \p\ cos + iO 

7T Q(-P°)Q(P 2 ) 1q / (-P°/2 + p°) 2 -(|fc| + |p|) 2 + i0 \ 
4 |P ||p| 0g l v (_po /2+:p O)2_ ( |^|_| jr | )2 + i0y ) 

7r 9(-p°)e(P 2 ) / p 2 -pV-|p||p°|+^o \ , , 

4 |P°||p| ^^-POpO + lpHPOI+iOy 1 j 

The distribution 8(Pq) inserted in eq. (|46|) does not affect the result in the rest frame of 
P, but will account below for the fact that I\ (p, q) vanishes for spacelike P when written in 
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covariant form @(P 2 ). Investigating the behavior of the expression |p||P°| under Lorentz 
transformations one finds that the covariant form of is given by V^V, where iV is 

the triangle function 

N = (Pp) 2 - P 2 p 2 = (pq) 2 

such that one arrives at 



„2 2 



p 



(47) 



h(p,q) 



ir&(-P°)Q(P 
4 y/N 

vr e(-P°)6(P 2 ) 




(48) 



It is left to the reader as a simple exercise to show that lightlike P do not contribute to 
I\(p,q) in a distributional sense. Calculating the remaining two integrals along the same 
lines as above and subtracting the resulting expressions with p and q interchanged, one 
obtains 



A 3 (p,9) 



vr J sgn(P°)9(P 2 



4(2vr) 



sgn(q°)e(q 2 ) 



N 



log 



p 2 — pq + \fN 



p z — pq 



'N 

sgn(p°)@(p 2 ) 



log 



pq + VN 



pq 



log 



N 
pq + 



g z — pq 



(49) 



Note that all imaginary parts of the integrals appearing in eq. (|40[) mutually cancel, since 
A3 is real due to its antisymmetry properties in configuration space described above. 



5. Calculation of A r 3 et 



In close analogy to eq. (|21j) . A r 3 et (p, q) is calculated in the following for p £ V + and q £ V . 
The result for arbitrary p, q can be obtained by analytic continuation, since A r 3 et (p, q) is 
the boundary value of an analytic function, regular in p G K 4 + iV + , q £ K, 4 + iV~ (see 
also HO]). The procedure is in two steps. First, the integral 



-00 



a 3p (p,<z) = - / dt^ML (so) 

—00 

is evaluated for p £ V + . This corresponds to the calculation of G(x? — x^)As(xi, X2, X3) in 
configuration space in the Lorentz system where p = (p°,0). Then A 3 et (p, q) is obtained 
by calculating 

+00 ~ +00 +00 

:retr„s_ * f A 3 , p (p, t?) 1 f M [ j_ A 3 (tp, S</) 



AS et (p,g) = — - / eft *£"2_^ = -— — / dt / ds 7 °\7 ^ -, (51) 

3 V ^' y; 2vr J 1-t + fO (2vr) 2 7 7 (1 - 1 + *0)(1 - a + tO) ' V ; 

—00 —00 —00 

corresponding to the multiplication of Q(x\ — x^)A^{x\, X2, x%) with 6(^2 — x®) in config- 
uration space. Note that the well-known splitting formula ^ 

+00 

.ret, ^ _ * /" M &a(tP,tq) 



A r , et (p,q) = — / tft ° v g ^ , (52) 
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is infrared divergent in the case of a totally massless vertex, while the double dispersion 
relation eq. (fSl"|) is not. 

For the actual calculations one may start by considering the first logarithm in eq. (|49jl 



+00 



h, P (p,q) 



ITT 



4(2vr) 



Sgnjfo - q°)e({tp - q) 2 ) 
at — — ^= log 



[i-t + mtWN 



tpq+ \t\VN 



tpq- \t\VN 



(53) 



The real part of the integral ()53j) can be written down immediately (P = p 



Re-[h, P (p,q)] 



tt 2 sgn(P°)e(P 2 ) 



4(2vr) 



N 



log 



pq + VN 



pq — y/N 



(54) 



The Lorentz system shall be chosen such that p — > (p°, 0). Since p G V + (and q arbitrary), 
one has to evaluate the Cauchy principal value 



Im[h, P (p, q)} 



ITT 



4(2vr) : 



log 



pq + VN 



+00 



dt 



sgn(tp° - q°)e(((tp°,6) - qf 
t(l-t) 



p<7 — y/N 

From the distribution, one obtains integral limits via 

t 2 p 2 - 2tp°q° + q 2 = 0^ t 1>2 = ^(q°± \A 2 - q 2 ) = ^(<7° ± 1*1), 



(55) 



(56) 



and taking into account that sgn(ti2p — q ) = sgn(±|<f|), one obtains from 

ta 



+ 



t(l-t) / t(l-t) 
((p°,0)- ? 2 



log 



log 



log 



(l-tl)(l-i 2 ) 



tlt 2 

(p - q) 2 



the first 'splitting' result 



h, P {p, q) 



ITT 



4(2vr) 3 ViV 



log 



pq + VN 



pq — V N 



log 



(p - qf 



vr 2 sgn(P°)e(P 2 



4(2vr) 



log 



(57) 
pq + VN 



As further step one may calculate 



+00 



1.P.9 



2tt 



h, P (p,tq) 
1-t + iQ' 



pq — y/ N 
(58) 



(59) 



The terms generated by the iO-term can be written down immediately 



ITT 



log 



pq + VN 



pq — y/N 



log 



(P - 



4(27r) 4 ViV 
The remaining non-trivial integrals are 



7T" 



4(27r) 4 ViV 



log 



pq+VN 



pq — y/N 



7T 



'1 



4(27r) 4 ViV 



log 



p<7 — \/N 



dt 



t(l-t) 



log 



(60) 



(61) 
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and 



iir 



1<2 



log 



pq + VN 



pq - VN 



dt 



■sgn(tP°)0(t 2 P 2 ^ 



t(l-t 



(62) 



Integral i\ can easily be solved by the help of the identity 

+R , I , 9 

hm / dt— — = — sgn(/3 - a), 

i?^oo J t — p I 



-R 



and one obtains 



7T ' 



n 



log 



pq + VN 



(63) 



(64) 



4(2tt) 4 \/]V " pq-y/N 

i.e. ii cancels the second term in eq. (|60[) . Integral %2 has a structure completely analogous 
to eq. (|55|) and is given by 



log 



pq + VN 



pq — \ N 



log 



4(27r) 3 VJV 

such the compact result (p £ V + , q £ V~) 

~ 2 1 'pq + y/N 



P 



(p - qf 



(65) 



I 



ITT 



log 



4(2^)4^ \pq-VN, 



log 



p 2 q 2 

(p - qf 



(66) 



is obtained. 

The integrals which appear in the calculation of the remaining parts of A3 (p, g), 
namely /2,p,ij and /3, p , g , can all be expressed by the help of the Spence function (or 
dilogarithm) Li(z), which is defined by 



Li(z) = J *log(l-* 



(67) 



where the integral contour in eq. Q67|) must be equivalent to a straight line from to z 
avoiding [1, +00) c II. This is due to the fact that all relevant integrals can be reduced 
directly to the form (a, b,c,d £ H) 



dt 



log |ai + 6| 1 j a(cx + d)\ 1 



ct + d 



--L 

c 



ad — be 



+ - log 

c 



ad — be 



log 



a{cx + d) 



ad — be 



where L(x) is the real Spence function, defined for real arguments x€H by 

L(x) = -(Li(x + iO) + Li(x - iO)). 
For the sake of brevity, the results are presented here 



(68) 



(69) 



h,p{p,q) 



ITT 



4(2vr) 



„ sgn(g°)e(g 2 ) 
dt- — — t= log 



(1 -t + iO)\t\VN 



t 2 p 2 -tpq+ \t\VN 



t 2 p 2 - tpq - \t\VN 
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ir 2 sgn((7°)0(g 2 



4(2vr) 



h, P ,q(p,q) 



N 



log 



p z — pq + \/N 



p* — pq 



ITT 



4(27r) 4 ViV 



2Li 



p — 



7T 



4(27r) 4 ViV 



log 



+ log 



— p 



2 (-pq + ^N" 



P 



log 2 



in 3 e(g°)e( g 2 )9(-p 2 

4(2^)3 
^ -2Li\ 
-pq — yN 



(70) 



p — pq 



N 



-pq + VN 



pq — p 2 + \//V 



+ log 



pg — q 



N 



pq — q 2 + VN 



log 



'pq + V^V N 



,pq 



(71) 



and /3, p ,q ist given by l3, p , q (p,q) = l2, P , q (q,p)- Note that the last three logarithmic terms 
in eq. (fTTj) are superfluous. They indeed appear in the actual calculations, however, they 
combine to zero, since 



' pq - p 2 - Vn\ ( pq - q 2 - Vn\ ( pq + VN^ 



pq — p 2 + V N J \pq — q 2 + \/N J \pq 



(72) 



It must also be pointed out that the Spence function obeys very many identities. E.g., the 
real Spence function satisfies 



L(l — x) = —L(x) + log |x| log |1 — x\ 



7T 

~6~ 



L(l/x) — — log 2 \x\ + log \x\ log |1 — x\ + 



-vr 2 /3 : x < 
+tt 2 /6 : x > 



-L(l - 1/x) + ^log 2 |x| + 



-vr 2 /2 : x < 
: x > 



(73) 



Looking for a compact expression for Ag e *, one finds 



27T 



2(2vr) 4 v / iV 



2Lf 



p — pq 



N 



+ log 



'gg_+ v^v' 



TV 



log 



2Li 



'p 2 - pq + v^/V^ 



P 



(P-9) 2 



(74) 



This expression can be transformed with the help of eq. ()73|) and the observation that the 
arguments of the Spence function can be written in many different ways 



p — pq 



N 



pq 



-pq - VN 
p 2 -pq+ VN 



-pq + V N 



1 



P 



pq + VN' 

■2 



P 



P<7 



iV 



(75) 
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Note that A^ et (p, q) is purely imaginary for p £ V + ,q £ V . This insight can also be 
obtained directly from the naive splitting formula eq. I|52j). where only the imaginary part 
is infrared divergent. There, one obtains for the real part of A^ et (p, q) (see also eq. (|71|l) 



2vH 



4(27r)ViV 



1q f pq-p 2 - VN \ t b f pq-q 2 - Vn \ ^ f pq + V^V ' 



0. (76) 



The retarded distribution A3 (p, (/) can be obtained for arbitrary p, q from eq. ()74|) by 
analytic continuation, since Ag e *(p, q) is the boundary value of the corresponding function 
analytic for p G H 4 + iV + and q G K 4 + i^ - . Furthermore, the full vertex function is 
finally obtained from eqns. 

y(p,g) = A- t (p, (? )-A3(p,g), (77) 

with 

A3 (p,g) = 




| e( g °)9( g 2 ) / p 2 - Pg + 

/]V 0g ^ p 2_^_^ A r * y ^ ~\q 2 -pq-y/N 

(78) 

The result for A 3 (p, g) is given here with the remark that A 3 (p, q) and A3(p, q) can indeed 
be calculated in parallel. A detailed calculation of A3 and A3 can also be found in pQ, 
chapter 3.8, where the vertex function appearing in QED with massive electrons is treated. 
The result for the massless case can be easily obtained by setting the electron mass equal 
to zero, since no infrared problems are present in A3 and A3. 

For p £ V + and q G V~ , A3 (p, q) vanishes and one can finally write down in a 'sloppy' 
style the result that can be found in the literature H3 H4j 



/ 



d 4 k 



[(p-k) 2 ][(q-k) 2 ][k 2 } 



ITT 2 

2yfN 



p 2 —pq—\fN\ . ( p 2 — pq + vN\ (pq + VN\i ( P 2 



2Li\ Z \-2Li\ c - +log - =7= log 



p 2 J \ p 2 J \pq - \fN J \(p — q) 2 

(79) 

Note that in the result given in 0, a factor of two is missing. Eq. ()79|) is also equivalent 
to the forms derived in |15l I16j . 
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